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PREFACE 


This  report  presents  results  of  a theoretical  study  which  evaluates 
the  response  of  hardened  underground  facilities  in  rock  to  dynamically 
applied  loads  produced  by  explosions.  This  research  was  conducted  by 
personnel  of  the  Phenomenology  and  Effects  Division  'PED)  of  the  Weapons 
Effects  Laboratory  (WEL),  U.  S.  Army  Engineer  Waterway.  Experiment  Sta- 
tion (WES),  during  the  period  January  1975-January  1976. 

The  primary  analytical  development  and  preparation  of  the  initial 
draft  report  were  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask 
J3kCAXSX311,  "Underground  Structures  Studies,"  under  the  guidance  of 
Dr.  Kent  Goering.  Final  report  preparation  and  additional  analytical 
work  including  investigations  of  backpacking  and  unloading  and  subse- 
quent reloading  were  sponsored  by  the  Office,  Chief  of  Engineers,  under 
Project  LA762719AT^0/Al/017,  "Stability  of  Deep  Underground  Structures 
in  Bock,"  which  was  monitored  by  Mr.  D.  S.  Reynolds. 

This  report  was  written  by  Messrs.  J.  L.  Drake  and  J.  R.  Britt, 

PED,  under  the  general  supervision  of  Mr.  L.  F.  Ingram,  Chief,  PED,  and 
Mr.  W.  J.  Flathau,  Chief,  WEL. 

COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE,  were  Directors  of  V?ES 
during  the  preparation  and  publication  of  this  report.  Mr.  F.  R.  Brown 
was  Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 

UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)  units  as  follows: 

Multiply To  Obtain 

mils  0. 0025^  centimetres 

inches  2.5^  centimetres 

pounds  (force)  per  6.89^757  kilopascals 

square  inch 

kilobars  100.00  megapascals 

inches  per  second  2.5^  centimetres  per  second 

per  second  per  second 

degrees  (angle)  0.017^5329  radians 
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A METHOD  FOR  DESIGNING  DEEP  UNDERGROUND 


STRUCTURES  SUBJECTED  TO  DYNAMIC  LOADS 


CHAPTER  1 
INTRODUCTION 


1 . 1 BACKGROUND 

Design  of  superhard  underground  facilities  to  resist  the  shock 
levels  induced  by  nuclear  weapons  must  take  full  advantage  of  the 
strength  of  the  surrounding  rock.  Survival  at  levels  in  excess  of 
2 kbars1  has  been  demonstrated  by  experiments  in  an  idealized  hard  rock 
medium.  In  principle,  additional  hardness  can  be  achieved  by  placing 
the  facility  deep  enough  to  attentuate  the  shock  to  levels  that  can  be 
resisted  by  structural  hardening  procedures. 

Static  methods  are  now  being  used  for  design  of  hardened  underground 
facilities.  Perhaps  the  best  analytical  procedure  is  that  outlined  by 
N.  M.  Newmark  (Reference  l)  wherein  rock  and  liner  systems  are  treated 
as  elastoplastic  materials  that  obey  a general  form  of  the  Mohr-Coulomb 
failure  criterion.  Incompressible  plastic  strains  and  axially  symmetric 
stresses  are  assumed.  This  analysis  proceeds  from  the  interior  of  the 
lining-medium  system  where  the  circumferential  strain  and  radial  stress 
are  known,  and  successively  works  the  solution  outward  by  matching 
stresses  and  displacements  across  liner  Junctions  and  finally  the  rock- 
liner  interface  to  obtain  the  free-field  conditions  in  the  rocV . Simi- 
larly, it  is  possible  to  work  from  known  conditions  on  any  interior  ele- 
ment inward  or  outward,  but  not  from  the  free-field  stress  situation. 

A.  J.  Hendron,  Jr.,  and  A.  K.  Aiyer  (Reference  2)  have  “xtended 
Newmark' s static  analysis  to  include  dilatancy  of  the  Mohr-Coulomb  mate- 
rial at  failure.  These  authors  report,  however,  that  the  increase  in 
volume  predicted  by  their  theory  is  too  large  compared  with  experimental 


A table  of  factors  for  converting  U.  S.  customary  units  of  measurement 
to  metric  (SI)  units  is  presented  on  page  3. 
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evidence  in  real  rocks.  They  claim  that  values  of  radial  displacement 
calculated  by  their  theory  are  a conservative  upper  bound.  On  the  other 
hand,  Newmark's  analysis,  which  neglects  any  volume  change  due  to  the 
plastic  strains  at  failure,  predicts  radial  displacements  that  are  too 
small  and  should  be  used  as  a lower  bound. 

Sophisticated  computer  codes  are  available  that  can  account  for  the 
complex  interactions  of  the  structural  elements  under  dynamic  loading. 

Gene  -al  stress-strain  relationships  that  closely  model  the  real  material 
behavior  in  both  the  elastic  and  inelastic  regions  can  be  used  in  theLe 
types  of  codes.  However,  these  procedures  are  still  too  complex  and 
time-consuming  to  be  useful  for  design  but  should  be  useful  for  final 
analysis  after  a preliminary  design  has  been  proposed. 

Treating  the  dynamic  load  as  an  equivalent  static  load  is  accept- 
able for  use  in  structural  design  methods.  In  the  case  of  normal  struc- 
tural elements  (columns,  beams,  etc.),  the  equivalent  static  loading  may 
be  determined  from  a dynamic  response  chart  for  various  ratios  of  the 
duration  of  loading  to  the  natural  period  of  the  structure.  For  rock 
tunnel  and  liner  systems  the  natural  period  is  not  easily  defined  with- 
out a dynamic  analysis  of  the  entire  system. 

1.2  OBJECTIVE 

The  goal  of  this  report  is  to  develop  solutions  for  a class  of  dy- 
namic problems  that  have  direct  bearing  on  design  considerations  for 
deep  underground  facilities  in  rock.  Specifically,  an  extension  of  the 
static  methods  outlined  by  Newnark  to  accommodate  dynamic  loads  is  sought. 

1.3  APPROACH 

This  report  describes  solutions  to  a dess  of  dynamic  elastoplastic 
problems  that  model  some  of  the  salient  features  of  the  response  of 
hardened  underground  facilities  in  rock.  The  theoretical  model  consists 
of  multilayered  concentric  cylinders  of  elastoplastic  materials  with 
time -dependent  loads  applied  to  the  exterior  boundary.  Each  element  in 
the  cross  section  is  assumed  to  be  incompressible  and  its  yield  governed 
by  a Mohr-Coulomb  failure  criterion.  The  number  of  elements  within  the 
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cross  section  is  not  limited.  Solutions  of  the  theoretical  model  are 
cast  in  the  general  form  normally  used  in  structural  dynamics: 

Mass  * Acceleration  = External  applied  load  - Internal  resistance 

The  resulting  equations  can  be  quickly  and  inexpensively  evaluated  on  a 
digital  computer. 

To  extend  the  range  of  validity  of  this  exact  theory,  a first-order 
correction  factor  was  developed  to  account  for  the  compressibility  of  the 
materials,  and  a simple  method  to  treat  backpacked  structures  was  intro- 
duced. The  theory  was  then  verified  by  comparing  calculated  values  with 
experimental  measurements  from  small-scale  static  and  explosively  driven 
tunnel  collapse  studies. 
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CHAPTER  2 


THEORY 


2.1  PROBLEM  FORMULATION 

The  rock  and  liner  system  is  modeled  as  multilayered  concentric 
cylinders  of  elastoplastic  materials  in  plane  strain  (see  Figure  2.1). 

An  arbitrary  time-dependent  pressure  is  applied  to  the  outer  boundary, 
ana  an  arbitrary  pressure,  also  time-dependent  if  desired,  is  applied  to 
the  internal  boundary.  It  is  assumed  that  the  materials  in  the  cross 
section  are  incompressible  in  both  the  elastic  and  the  plastic  states. 
(As  discussed  in  Section  2.3,  a first-order  compressibility  correction 
is  easily  incorporated.)  With  these  assumptions,  the  effects  of  mate- 
rial compaction  and  stress  wave  interaction  within  the  cross  section  are 
neglected,  thereby  greatly  simplifying  the  analysis.  The  number  of  ma- 
terial layers  in  this  analysis  is  unlimited. 

2.1.1  Field  Equations.  Equations  requiring  equilibrium  and  conser- 
vation of  mass  are  basic  to  an}  analytical  study  of  structural  response. 
Solutions  to  these  equations  are  valid  in  both  the  elastic  and  plastic 
states  of  the  material. 

The  equation  of  motion  in  polar  coordinates  is  as  follows:^ 


i 1 1 1 /-2  : 

do  a - o„  / 3 u 

r . r 6 _ ( i 

3r  r " Pi  Vat2 


!i  = 1,2,. . . ,m) 


where 


o = radial  stress  component 
r 

r = radial  coordinate 
Og  - tangential  stress  component 
= mass  density 

ui  = radial  displacement  component 
t - time 


For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  A). 


UNIFORM  EXTERIOR  PRESSURE  pit) 


t 


LEGEND 

<2  = INTERIOR  BOUNDARIES  OF  THE  FIRST  AND 
SECOND  MATERIAL  LAYERS.  RESPECTIVELY 

HI  - NUMBER  OF  MATERIALS  IN  THE  CROSS 
SECTION 

P - DENSITY 

G = SHEAR  MODULUS 

q - COHESION 

(f,  - ANGLE  OF  INTERNAL  FRICTION 

Figure  2.1  Geometry  of  rock-liner  system. 


i = subscript  or  superscript  referring  to  a component  material  in 
the  cross  section 

m = number  of  materials  in  the  cross  section 

This  equation  can  be  solved  when  relationships  between  stresses  and  dis- 
placements (constitutive  laws)  are  given.  In  the  plastic  state  an  addi- 
tional equation  given  by  the  yield  criterion  relating  the  radial  and 
tangential  stresses  is  necessary  for  a complete  solution. 

Incompressibility  requires  the  conservation  of  mass  to  be 


3U1  u1 
— L + -L  = o 
3r  r 


(i  = 1,2 :n) 


(2.2) 


The  general  solution  to  Equation  2.2  is 


„i  ,thl 

r r 


(2.3) 


where  U^t)  is  a general  function  of  time  tc  be  determined  from  the 
boundary  conditions.  It  is  evident  from  Equation  2.3  that  if  the  dis- 
placement field  is  continuous  between  component  materials,  U*(t)  must 
be  the  same  for  all  materials.  Thus,  the  superscript  i is  dropped 
such  that 


U!(t)  = u2(t)  =...=  U(t) 


(?.M 


2.1.2  Material  Models.  Stress-strain  relationships  for  an  incom- 
pressible elastic  medium  in  plane  strain  can  be  written  as 


o1  = S1  + 2G . e 1 
r l r 


°9  = s‘  * xice 


1 - e1 
0 = 0 


(2.5) 


wh  “re 


s1 


u . 

l 


mean  normal  stress 
shear  modulus 
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e = radial  strain  component 

= tangential  strain  component 
= axial  stress  component 
and  the  mean  normal  stress  is  given  by 


-i  _ 1 (i 

3 ' 3 [°r 


+ o ♦ a 

0 z i 


In  polar  coordinates,  the  strains  are  defined  as 

i 


3u 

l _ r 

Er  = 3r 


U(t) 

2 


i _ ^r  _ U(t) 

E0  ' r 2 

r 


(2.6) 


where  the  radial  displacement  given  by  Equation  2.3  is  substituted  into 
the  strain  definitions.  Because  U(t)  is  independent  of  the  state  of 
the  material,  it  is  evident  that  Equations  2.6  are  valid  in  both  the 
elastic  and  plastic  regions. 

Stress-strain  relations  in  the  plastic  region  are  implied  and  de- 
fined by  the  condition  of  incompressibility  and  from  the  form  of  the 
yield  law.  A Mohr-Coulomb  yield  criterion  assumed  for  all  the  component 
materials  is  given  as 


(2.7) 


wiit‘  re 


N. 


1 + sin 
1 - sin 


and  q.  is  the  cohesion  and  is  the  angle  of  internal  friction. 

The  assumption  of  incompressibility  is  not  strictly  compatible  with  this 
yield  c»  iterior.  except  for  <t.  = 0 (von  Mises  yield  criteria)  unless  a 
nonassociative  flow  rule  is  assumed.  Mohr-Coulomb  materials  dilate 
while  undergoing  failure;  however,  neglecting  this  effect  is  not 
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expected  to  affect  the  results  greatly,  and  simplifies  the  solution. 

2.1.3.  Boundary  Conditions.  Stresses  and  displacements  are  assumed 
to  be  continuous  across  the  various  internal  boundaries  between  mate- 
rials. Further,  dynamic  pressures  are  applied  at  the  exterior  and  in- 
terior boundaries.  The  external  pressure  was  assumed  to  be  much  greater 
than  the  interior  pressure,  causing  the  deformation  to  be  inward.  These 
boundary  conditions  can  be  expressed  as 


o 


1 


o 

u 

o 


r 

i 

r 

i 

r 

m 

r 


-Po(t) 


-p(t) 


at  r = r^ 
at  r = ri.i 


at  r 


rm+l 


(i  = 1,2,.. . ,m-l) 


(2.8) 


where 

p (t)  = internal  pressure  (compressive) 

i\  = interior  boundary  of  the  i material  layer 
p(t)  = external  pressure-time  history  (compressive) 
all  of  which  are  depicted  in  Figure  2.1. 

Additional  internal  bou”  ..ary  conditions  may  occur  if  any  of  the  com- 
ponent materials  are  partially  elastic  and  partially  plastic.  Across 
this  elastic-plastic  boundary 


(0r"°e)\  =(°r-°e)\ 

elastic  plastic 


and  the  normal  stresses  o are  continuous, 

r 


2.2  PROBLEM  SOLUTIONS 

Solu* ions  to  the  equation  of  motion  (Equation  2.1)  are  possible 
from  the  displacement  solution  (Equation  2.3)  and  either  the  stress- 
strain  relations  in  the  elastic  region  (Equation  2.5)  or  the  yield 
criterion  (Equation  2.7)  in  the  plastic  region. 

2.2.1  Elastic  Region.  In  the  elastic  region,  the  following 
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equation  is  obtained  from  the  elastic  stress-strain  relations  and  from 
Equations  2.6: 


■;  - °e  * --i|¥ 

_ r 


Substituting  this  expression  into  the  equation  of  motion  (Equation  2.1) 


gives 


with  the  solution 


i , /r  1 d“"U(t) 

a = p . lnj — — 

r 1 Vi  dt2 


2G . /r 


U(t)  + C. 


(2.9) 


where  is  a constant  independent  of  r to  be  determined  from  bound- 

ary conditions.  From  this  result  and  Equations  2.5  and  2.6  the  mean 


normal  stress  S is  given  by 


/2G.\/r.\2 

S1  = o1  + U(t) 

r I J\ r J 


(2.10) 


2.2.2  Flastic  Region.  An  expression  for  the  stress  difference 

in  the  plastic  region  is  given  by  the  yield  condition  (Equa- 
tion 2.7).  Substituting  this  relationship  into  the  equation  of  motion 
(Equation  2.1)  yields 


with  the  solution 


‘'j-#' 
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Note  that 


n.-*-  o 
1 


i&f-Hy 
■ - 


where  n.  = N.  - 1 . Thus,  in  the  limit  as  N.  -*  1 or 
l l ’ i 


-*•  0 , the  ex- 


pressions for  the  elastic  stress  distribution  (Equation  2.9)  and  plastic 
distributions  (Equation  2.1l)  bear  similar  terms  for  the  inertial  con- 
tribution and  for  the  arbitrary  constant  C\  . 

2.2.3  Generalized  Stress  Distribution.  Since  each  material  is  al- 
lowed to  be  elastic  or  plastic  or  even  partially  yielded,  the  number  of 

stress  combinations  required  to  solve  the  boundary  conditions  for  the 

2 

constant  C.  is  m for  m materials  in  the  cross  section.  This  sit- 
l 

uation  can  be  avoided  by  writing  a more  general  stress  distribution  which 
can  be  used  in  both  the  elastic  and  plastic  regions.  This  generalized 
stress  distribution  is  expressed  in  operator  form  as 

\i-  W )}" 


- 1 ♦ c . — 


(2.12) 


where 


2 s^n 

n.  = N.  - 1 = t : — - 

ii  1 - sm  d 


and  again  noting 


- 1 = In 


In  the  elastic  state  it  is  necessary  to  define  n^  = q^  =0  for  Equa- 
tion 2.12  to  be  valid.  The  plastic  state  requires  = 0 to  obtain  the 
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correct  stresses.  Using  Equation  2.12  in  the  boundary  conditions  will 
yield  a generalized  form  for  the  constant  which  is  valid  in  either 

the  elastic  or  plastic  state  and  which  correctly  accounts  for  the  states 
of  all  m elements. 

2.2.14  Stress  Solution.  Substituting  Equation  2.12  into  the  boundary 
conditions  expressed  by  Equations  2.8  gives  the  following  recurrence 
formulas  for  C.  : 

l 


C1  * -Po(t>  * — “ 

rl 


c .fipa 

i+1  n.  I r. 

i \ i 


.2..  20. /r.  \ 2G.  U 

dJJ  _ — LI — L_ | u + _1+L_ 
,2  2 [r.J  U 2 

dt  r.  \ 1+1/  r... 

l \ / l+l 


* n. 

2q . /r. . ,\  1 

— ) - 

n.  y l r, 

l \ l J 


1 + c, 


(2.13) 


Completing  this  solution  and  continuing  the  recurrence  out  to  r = rm+1 
gives  the  following  formula: 


■/  vn.  "|  , .n.  } 0 

!i  (t_\  * . , . /s_V*  (A 

"i[\ri/  J VV  -‘jdt2 

icG.  T/  \"i*2  1 / \"i  ) 

•J/p)  - 'J  * fc)  "."I" 

■ ' ']  * (s?)  Vi..  • 


i 


where 


i 


where 


J=1  J |/  J Jk«j+l'  ' 

, yi/^^rrrr  [fr±4  J - ill  fT 

1 bhr 3 IW  JLvArJ 


i 30.  L...X  l J. 


- nfe) 

J>1  J*1  ' J ' k=J*l  ' ' 


Poi  * "o'1 11 


np?) 


>(4 


PQ(t)  - pQ(t)  - (KU  - Y) 


(2.15) 


(2.l6) 


M = M * effective  mass 
m 

p * p s effective  interior  pressure 
o om 

K ■ K * effective  elastic  stiffness 
m 

Y * Y » effective  rigid  plastic  collapse  pressure  of  the  m 
materials 


The  symbol  ) denotes  a sum  of  terms  with  J * 1 to  J * i , and 


1 l denotes  a product  of  terms  with  k = J + 1 to  k = i . 
k*J+l 

Equations  2.1U-2.16  along  with  the  yield  condition  (Equation  2.7) 
completely  describe  the  stresses  and  motions  of  the  rock-liner  system 
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J. 


under  dynamic  loading.  The  differential  Equation  2.l6  is  of  the  same 
form  as  that  used  in  structural  dynamics: 

Mass  x Acceleration  = External  applied 
load  - Internal  resistance 


In  this  case,  the  internal  load-resistance  function  is  R = KU  - Y . 

Note  that  initially  Y = 0 and  the  resistance  is  due  to  the  elastic 
stiffness  K . As  the  external  applied  load  increases,  elements  begin  to 
fail,  effectively  increasing  Y as  K decreases  until  the  applied  load 
is  such  that  the  entire  section  fails,  the  full  value  of  Y is  reached, 
and  K = 0 . 

2.2.5  Elastic-Plastic  Boundaries.  Solutions  expressed  by  Equa- 
tions 2.11+  and  2.16  are  valid  for  all  load  conditions  in  both  the  elas- 
tic and  plastic  ranges  as  long  as  the  deformations  are  inward;  however, 
these  equations  alone  are  not  sufficient  to  describe  the  stresses  and 
motions  of  the  rock-liner  system.  Yield  conditions  must  be  used  along 
with  these  equations  to  determine  the  elastic-plastic  state  of  the 
system. 

When  the  deformation  is  inward,  yielding  in  each  of  the  component 
materials  begins  at  the  interior  boundary  and  progresses  outward  until 
the  entire  element  is  fully  plastic.  This  behavior  can  be  used  to 
simplify  the  analysis  as  follows.  Divide  each  component  material  into 
two  parts,  the  inner  portion  with  only  plastic  properties  and  the  outer 
shell  with  elastic  properties.  With  this  definition  the  m materials 
are  divided  into  two  m layers  where 


n 

q 

G 


i 

i 

i 


0 

0 


for  even  values  of  i 
for  odd  values  of  i 


(2.17) 


and  r^+^  for  odd  values  of  i is  the  elastic-plastic  interface.  A 
material  is  fully  elastic  when  r^+^  * for  odd  values  of  i and 
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a material  is  fully  plastic  when  for  odd  values  of  i . 
Using  these  definitions,  the  elastic-plastic  boundary,  r^+^  for  °dd 
values  of  i , can  be  calculated  from  the  yield  conditions. 

Continuity  of  stresses  across  the  elastic-plastic  interface 
requires  that 


/ i+1  i+l\  / i i\ 
\°r  - °e  ) = (ar  - °ej 


at  r = r^+^  for  odd  values  of  i 


Using  Equations  2.5  and  2.6  in  the  elastic  range  and  Equation  2.7  in  the 
plastic  region  gives 

2 


- -jr1  U * -Vr  * 2qiVnl  ♦ 1 at  r ■ r 


i+1 


for  odd  values  of  i 


Substituting  in  this  result  the  solution  for  radial  stresses  (Equa- 

I 

tion  2.lU)  yields  the  expression  for  r ,+^/r^  and  i - odd 
« n, 


i r 


(p.  + n.M.  ,)f4' 
11  1 1 1-1'\dt2; 


+ n.^K.  U - Y.  . - P . , ) - 2q.  <Jn.  + 1 
i\  1-1  1-1  01-1/  1 \ 1 


+ (2>l8) 
\dt  / 


For  n^  = 0 this  equation  reduces  to  the  simple  form 


Pf)  • 

\ i / Vi 


(2.19) 


Note  that  for  the  Mohr-Coulomb  failure  criterion,  inertial  mass  * 
acceleration  terms  appear  in  the  condition  for  dynamic  failure.  Until 


the  time  at  which  the  maximum  inward  particle  velocity  is  reached,  when 


2 2 

d U/dt  0 , the  inertial  forces  tend  to  strengthen  the  section  causing 

o 


r.  ,/r.  to  be  less  than  its  static  value.  However,  when  d“~U/dt  > 0 , 
i+l  l 
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these  forces  cause  the  plastic  boundary  to  grow  outward  more  rapidly 
than  under  static  load. 


2.3  A FIRST-ORDER  COMPRESSIBILITY  CORRECTION 

For  static  loading  and  elastic  response  the  ratio  of  the  compres- 
sible to  the  incompressible  radial  displacements  is  (Reference  3) 


u ? 

^-  = 1 + (1  - 2v)r 


r2P2  - rlPl 
rlr2(p2  - Pl> 


where  v is  Poisson's  ratio  of  the  material,  r and  p^  are  the  in- 
side radius  and  applied  loa  . respectively,  and  r^  and  p^  are  the 
outside  radius  and  applied  j.cad,  respectively.  For  most  tunnel  liners 
being  considered,  p = 0 (atmospheric  pressure)  so  that  P-,/p0  is 
very  small.  Expanding  the  expression  for  u^u^  evaluated  at  r * ^ 
in  powers  of  p^Pg  gives 


1 + (1  - 2v)  <1  + 


, .2]  » , .n 


For  p^/p^  = ® this  expression  becomes 


u (r.) 

= 2(1  " v) 


(2.20) 


When  p^  is  zero,  this  ratio  is  exact  for  the  static  elastic  case. 

Experience  indicates  that  this  is  also  a reasonably  accurate,  first-order 

correction  factor  for  the  dynamic  displacements  of  an  elastic-plastic 

rock-tunnel  liner  system.  In  this  application  an  average  value  of  the 

Poisson's  ratios  of  the  component  materials  should  be  used.  Typically, 

v for  common  rock  and  liner  materials  ranges  from  about  0.2  to  0.3,  and 

uc(r^)/u^(r^)  ranges  from  1.6  to  1.1*.  Average  values  of  v * C.25  and 

u (r, )/u. (r, ) * 1.5  are  adequate  for  most  calculations, 
cl  i 1 
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2.14  A SUMMARY  OF  THE  CALCULATIONAL  METHOD 
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Equations  2.l8  or  2.19  together  with  the  equation  of  motion,  2.l6, 
can  be  solved  for  U(t)  by  a numerical  integration  such  as  the  Runge- 
Kutta  method.  Then  using  Equation  2.6  and  the  compressibility  correc- 
tion, Equation  2.20,  the  radial  displacement  can  be  determined.  The 
radial  stress  component  is  given  by  Equations  2 . 1^4  and  2.15;  then  with 
the  mean  normal  stress  S1  given  in  Equation  2.10,  the  other  stress 
components  are  given  in  Equations  2.5  and  2.7.  Thus,  these  relations 
provide  a complete  description  of  the  stresses  and  motions  of  the  rock- 
liner  system  in  a form  suitable  for  fast  computer  calculations. 

2.5  TREATMENT  OF  BACKPACKED  STRUCTURES 

Backpacking  is  commonly  a crushable  material,  such  as  foamed  con- 
crete, placed  between  the  rock  and  the  tunnel  liner  to  absorb  the  large 
strains  suffered  by  the  rock  and,  consequently,  to  isolate  the  liner 
from  excessive  damage.  The  simplest  treatment  of  backpacking  is  the 
hydrostatic  model,  which  neglects  the  shear  strength  of  the  material. 

Thus  the  pressures  at  the  inside  and  outside  of  the  backpacking  are  the 
same.  Experimental  pressure-volumetric  strain  curves  car.  be  used  in 
this  model  to  calculate  the  pressure  on  the  rock  and  liner.  A simple 
iterative  process  can  be  set  up  which  links  the  radial  displacement  of 
the  liner  to  the  displacement  of  the  rock  as  a function  of  volume  changes 
in  the  backpacking.  Since  the  present  analysis  allows  a time-varying 
interior  pressure  P()(t)  , the  calculational  method  is  easily  adapted 
for  computations  using  this  hydrostatic  model  of  backpacking.  The  com- 
puter code  presently  in  use  at  the  U.  U.  Army  Engineer  Waterways  Experi- 
ment Station  includes  this  option. 

For  free-field  stresses  in  the  rock  which  are  essentially  axisym- 
metric  (uniform  radial  load),  the  hydrostatic  model  provides  an  accurate 
description  of  the  backpacking  collapse.  The  effect  of  shear  strength 
of  the  backpacking  on  rock-liner  displacements  is  negligibly  small  in 
this  case. 

For  the  biaxiul  free-field  stresses  in  which  the  difference  between 
the  vertical  and  horizontal  normal  stress  components  is  large  compared 
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with  the  collapse  strength  of  the  backpacking,  the  hydrostatic  model  of 
collapse  of  the  backpacking  is  no  longer  valid  for  relatively  large  rock 
deformations.  Since  the  material  cannot  flow  to  equalize  the  stresses, 
ovaling  of  the  backpacking  results,  transmitting  highly  nonuniform  radial 
loads  to  the  liner  which  may  collapse  prematurely. 

2.6  DISCUSSION  OF  UNLOADING 

The  present  analysis  is  limited  to  monotonically  increasing  dis- 
placements. In  the  case  of  small  plastic  strains,  unload  and  a subse- 
quent reload  may  take  place  elastically  which  can  be  treated  simply 
within  the  present  analysis.  For  von  Mises  materials  in  static  loading, 
Prager  and  Hodge  (Reference  3)  give  equations  for  both  elastic  and  plas- 
tic unload.  The  plastic  unload  of  a Mohr-Coulomb  material  may  be  treated 
in  a similar  manner;  but  the  analysis  is  more  complex  than  the  von  Mises 
case,  as  the  authors  of  the  present  report  have  found  in  a thorough  study 
of  this  process.  In  the  present  study  it  was  discovered  that  the  Mohr- 
Coulom''  plastic  unloading  often  involves  further  yielding  of  the  material 
even  as  the  inward  displacement  decreases.  The  manner  and  extent  of  this 
unload  yielding  depend  primarily  on  (l)  the  amount  of  plastic  deformation 
in  the  material  at  the  end  of  the  load  phase  and  (2)  the  confining  pres- 
sure at  the  interior  of  the  layer  of  material.  The  resulting  equations 
are  mathematically  too  cumbersome  to  be  included  in  the  analysis  of  this 
report.  Th"  experience  gained  in  this  study  indicates  that  elastic  un- 
load is  often  a good  first  approximation  for  the  purposes  of  this  report, 
but  an  accurate  theoretical  description  of  repeated  inelastic  loading 
and  unloading  of  Mohr-Coulomb  materials  most  be  more  rigorous  and  will 
require  further  theoretical  work. 
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CHAPTER  3 


EVALUATION  OF  THE  THEORY 

3.1  COMPARISON  OF  CALCULATIONS  WITH  EXPERIMENTAL  EVIDENCE 

Calculations  were  made  for  small-scale  tests  conducted  by  the 
Stanford  Research  Institute  (SRI)  (Reference  k)  to  evaluate  the  theory 
of  this  report.  In  these  experiments  explosives  were  used  in  a spe- 
cially designed  device  to  produce  axially  symmetric  loads.  Displace- 
ments as  a function  of  time  were  measured  using  the  "light  ring"  tech- 
nique which  consists  of  photographing  the  motion  of  the  inner  pipe  wall 
by  means  of  a bright  ring  of  light  transmitted  through  a transparent 
flexible  plastic  tube  glued  to  the  pipe  wall.  The  targets  considered 
for  calculation  were  3-inch-OD  tubes  of  "super  lean  grout"  (SLG)  lined 
with  a 10-mil-thick  pipe  of  Type  3^7  stainless  steel . The  material 
properties  used  in  the  calculations  are  given  in  Table  3.1. 

Figure  3.1  shows  data  from  two  shots  with  80  grams  of  explosive. 

The  estimated  free-field  pressure  pulse  obtained  from  calibration  tests 
is  also  shown  in  the  figure.  Calculated  curves  are  for  Poisson's  ratio, 
v = 0.25  , 0.38  , and  0.5  , where  O.38  was  the  value  of  v calculated 
from  ultrasonic  velocity  measurements  in  SLG  given  in  Reference  5. 

There  is  good  agreement  between  the  theory  and  experiment.  Figure  3.2 
presents  measurements  from  three  shots  with  estimated  free-field  pres- 
sure pulses  as  shown.  Here  the  calculated  values  are  close  to  the  ex- 
perimental points,  but  the  agreement  is  not  as  good  as  in  Figure  3.1. 

An  additional  check  of  the  theory  can  be  obtained  by  noting  that 
the  internal  resistance  R = KU  - Y is  the  static  stress-strain  curve. 
Thus  data  from  static  or  quasi-static  tests  can  be  used  to  evaluate  this 
important  part  of  the  dynamic  equation  of  notion  (Equation  2.l6).  Fig- 
ure 3.3  shows  static  calculations  compared  v h measurements  obtained 
by  SRI.^  In  this  test  a hydrostatic  load  was  oplied  to  a sample  of 

^ Letter  of  10  Jun  1975  from  T.  C.  Kennedy,  Stanford  Research  Institute, 
Menlo  Park,  Calif.,  to  J.  Drake,  Weapons  Effects  Laboratory,  U.  S. 

Army  Engineer  Waterways  Experiment  Station,  CE,  Vicksburg,  Miss. 
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ime  for  three  press  ire  pulses  (from 


Figure  3.3  Computed  and  measured  static  tunnel 
closure  versus  pressure. 

WES  6B  rock  simulant  containing  a 0.625-inch-ID  tunnel  lined  with  a 
0 . 025-inch-thi ck  mild  steel  tube.  The  material  properties  used  in  the 
calculation  are  given  in  Table  3.1.  The  displacement  or  tunnel  closure 
is  expressed  as  a percentage  of  the  inside  diameter  D of  the  tunnel. 

The  calculated  unload  curves  assumed  elastic  behavior.  Again,  there  is 
good  agreement  between  theory  and  experiment. 

3.2  CONCLUSIONS 

An  analysis  has  teen  presented  which  extends  the  static  methods 
outlined  by  Newmark  (Reference  l)  to  accommodate  dynamic  loads.  Compari- 
sons of  calculations  with  experimental  data  have  demonstrated  that  the 
analysis  will  predict  satisfactorily  the  time-displacement  curve  in  uni- 
form radial  load  for  materials  which  can  bs  modeled  as  elastic-plastic 
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media  having  a Mohr-Coulomb  or  von  Mises  yield  criterion.  The  analysis 
given  in  this  report  is  in  a form  which  allows  quick,  inexpensive  pre- 
liminary design  calculations  for  deep  underground  structures. 

Users  of  the  theory  must  be  warned,  however,  that  the  analysis 
will  not  predict  accurately  displacements  for  biaxial  loads  where  the 
vertical  and  horizontal  normal  stress  components  are  significantly  dif- 
ferent or  where  a large  amount  of  ovaling  of  the  tunnel  liner  may  be 
expected.  In  this  case  the  theory  will  provide  a conservative  lower 
bound  estimate  of  the  deformation.  At  present,  a design  method  for 
biaxial  loading  is  not  available  in  a relatively  simple  form  comparable 
with  the  analysis  procedure  outlined  in  this  report.  Work  in  this  area 
is  needed. 

The  present  analysis  is  limited  to  monctonically  increasing  dis- 
placements. In  the  case  of  small  plastic  strains,  unload  and  a subse- 
quer * reload  may  take  place  elastically.  This  process  can  be  treated 
‘-imply  within  the  present  analysis. 
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TABLE  3.1  MATERIAL  PROPERTIES  USED  IN  THE  CALCULATIONS 


Material 

Shear 

Modulus  G 
6 

10  psi 

Poisson' 
Ratio  v 

Super  lean 
grout 

0.058 

0.38 

3bj  stainless 
steel 

11.2 

0.3 

WES  6b  rock 
simulant 

0.1+8 

0.25 

Mild  £teel 

12 

0.3 

Cohesion  q 
10^  psi 

Angle  of 
Internal 
Friction  <J> 
degrees 

Density  p 

-h 

10  psi 

in. /sec 

0.016 

0.7 

1.65 

20 

0 

7.5 

1.7 

30 

— 

20 

0 
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APPENDIX  A:  NOTATION 


G. 
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K 

m 

M 


n. 

1 

N. 
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p(t) 

P0(t) 

Po(t) 


PI,P2 


*i 

r 

r. 

l 

rl’r2 

R 

ni 


u /u. 
c 1 


Vue 

u 

Y 


i i 
cr’ee 


°r ,O0’°2 


An  integration  constant  in  the  radial  stress  solution 
Shear  modulus 

A subscript  or  superscript  referring  to  a component  material 
in  the  cross  section 

Effective  elastic  stiffness 

Number  of  materials  in  the  cross  section 

Effective  mass 

N.  - 1 

l 

(l  + sin  ) /l  - sin  4k) 

Exterior  pressure-time  history  (compressive) 

Interior  pressure  (compressive) 

Effective  interior  pressure  in  equation  of  motion 

Inside  and  outside  applied  load,  respectively 
Cohesion 

Radial  coordinate 

Interior  boundary  of  the  i material  layer 
Inside  and  outside  radii,  respectively 
Internal  resistance 
Mean  normal  stress 
Time 

Ratio  of  the  compressible  to  the  incompressible  radial 
displacements 

Radial  and  tangential  displacement  components,  respectively 
Displacement  potential 

Effective  rigid  plastic  collapse  pressure  of  the  m 
materials 

Radial  and  tangential  strain  components,  respectively 

Poisson's  ratio 

Density 

Radial,  tangential,  and  axial  stress  components,  respectively 
Angle  of  internal  friction 
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90  Church  Street,  New  York,  N.  Y.  10007 
ATTN:  NADEN/Chief,  Engineering  Division 

Division  Engineer 

U.  S.  nrmy  Engineer  Division,  North  Central 
536  South  Clark  Street,  Chicago,  111.  60605 


1 
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Address 


District  and  Division  Offices  (Continued) 
Division  Engineer 

U.  S.  Army  Engineer  Division,  North  Pacific 
210  Custom  House,  Portland,  Oreg.  97209 

Division  Engineer 

U.  S.  Army  Engineer  Division,  Ohio  River 
P.  0.  Box  1159,  Cincinnati,  Ohio  1+5201 

Division  Engineer 

U.  S.  Army  Engineer  Division,  Ohio  River 
P.  0.  Box  1159*  Cincinnati,  Ohio  1+5201 
ATTN:  Division  Laboratory,  0RDED-FL 

Division  Engineer 

U.  S.  Army  Engineer  Division,  South  Atlantic 

510  Title  Building,  30  Pryor  Street,  SW.  Atlanta,  Ga.  30303 

Division  Engineer 

U.  S.  Army  Engineer  Division,  South  Pacific 

630  Sansome  Street,  Room  1216,  San  Francisco,  Calif.  91+111 

Division  Engineer 

U.  S.  Army  Engineer  Division,  Southwestern 

Main  Tower  Building,  1200  Main  Street,  Dallas,  Tex.  75202 

Army 

Defense  Documentation  Center 

Cameron  Station,  Alexandria,  Va.  ?23ll+ 

ATTN:  TC/Mr.  Myer  B.  Kahn 

Director 

U.  S.  Army  Coastal  Engineering  Research  Center 
Kingman  Building,  Fort  Belvoir,  Va.  22060 
ATTN:  Librarian 

Commander /Director 

U.  S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
P.  0.  Box  282,  Hanover,  New  Hampshire  03755 
ATTN:  Library 

Director 

U.  S.  Army  Construction  Engineering  Research  Laboratory 
P.  0.  Box  1+005,  Champaign,  111.  61820 

ATTN:  Library 


No.  of 
Copies 
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1 


1 

1 


1 
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Copies 


.ontmuec 


Commander 

U.  S.  Army  Mobility  Equipment  R&D  Center 
Fort  Belvoir,  Va.  22060 

ATTN:  Technical  Documents  Center,  Building  315 

Commandant 

U.  S.  Army  Engineer  Scnool 
Fort  Belvoir,  Va.  22060 
ATTN:  Department  of  Engineering 

Secretary 

Commander 

U.  S.  Army  Mobility  Equipment  i-.&D  Center 
Fort  Belvoir,  Va.  22060 
ATTN:  Technical  Library 

Commander 

U.  S.  Army  Missile  Command 
Redstone  Scientific  Information  Center 
Redstone  Arsenal,  Ala.  35809 
ATTN:  Chief,  Documents 

Commander 

d.  S.  Army  Forces  Command 
Fort  McPherson,  Ga.  30330 
ATTN : AFEN-FEB 

AFEN  (Mil  Engr  Div) 

138th  Engineer  Group  (Construction) 

Department  of  the  Army 
Fort  Riley,  Kans.  66hk2 
ATTN:  ALBFEN-0 


Commandant 

U.  S.  Army  Command  and  General  Staff  College 
Fort  Leavenworth,  Kans.  66027 
ATTN:  ATSW-RI-L 


Commander 

Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adelphi , Md.  20783 

ATTN:  AMXDO-NP,  Mr.  Louis  Beiliveau/Mr . J.  H.  Gwaltney 
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Copies 


Army  (Continued) 

Director 

U.  S.  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Md.  21005 
ATTN:  AMXBR-X/Mr.  J.  J.  Meszaros 

Commander 

U.  S.  Army  Electronics  Command 
Fort  Monmouth,  N.  J.  07703 
ATTN : AMSEL-GG-TD 

Commanding  Officer 
Picatinny  Arsenal 
Dover,  N.  J.  07801 
ATTN:  ORDBR-TK 

Chief,  Plastics  Technical  Evaluation  Center 
SARPA-TS-S/No.  59 
P.  Angelloti 
SMUPA-ND-S/Mr.  E.  Zimpo 

Superintendent 
U.  S.  Military  Academy 
West  Point,  N.  Y.  10996 
ATTN : Library 

U.  S.  Military  Academy 
Department  of  Engineering 
West  Point,  N.  Y.  10996 
ATTN : MADN-F 

President 

U.  S.  Army  Air  Defense  Board 
Fort  Bliss,  Tex.  79916 

Commandant 

U.  S.  Army  Air  Defense  School 
Fort  Bliss,  Tex.  79916 

Commander 

U.  S.  Army  Nuclear  Agency 
Fort  Bliss,  Tex.  79916 
ATTN : ATCANA-W 

Director 

Eustis  Directorate,  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory 
Fort  Eustis,  Ya.  2360L 
ATTN:  SAVDL-EU-M0S  (Room  209) 
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Army  (Continued) 


Commander 

U.  S.  Army  Training  and  Doctrine  Command 
Fort  Monroe,  Va.  23351 

ATTN : ATEN-ME  1 

Commandant 

Armed  Forces  Staff  College 
Norfolk,  Va.  23511 

ATTN : Library  1 

Commander 

U.  S.  Army  Materiel  Development  and 
Readiness  Commani  (DARCOM) 

5001  Eisenhower  Avenue 
Alexandria,  Va.  22333 

ATTN : DACRD-WN  2 

Office,  Chief  of  Research,  Development, 
and  Acquisition  ( DAMA-CSS-D ) 

Department  of  the  Army 
Washington,  D.  Z 20310 
ATTN:  Mr.  Paul  F.  Carlton 

Technical  Library 

Navy 

Commanding  Officer  and  Director  1 

U.  S.  Naval  Electronics  Laboratory 
San  Diego,  Calif.  92152 

Commander  1 

Naval  Weapons  Center 
China  Lake,  Calif.  93555 

Superindent 

U.  S.  Naval  Postgraduate  School 
Monterey,  Calif.  939^0 

ATTN : Lib'-ary,  Code  212 1 1 

Commander 

Naval  Weapons  Center 
China  Lake,  Calif  93555 
ATTN:  Peggy  Davis,  Code  1*503 
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Address 


No.  of 
Copies 


Navy  (Continued) 


Officer  in  Charge 
Civil  Engineering  Laboratory 
Naval  Construction  Battalion  Center 
Port  Hueneme,  Calif.  930b3 

ATTN:  Code  L31  1 

Commanding  Officer  1 

Nuclear  Weapons  Training  Group,  Pacific 
U.  S.  Naval  Station,  North  Island 
San  Diego,  Calif.  92135 

Commanding  Officer  and  Director 

Naval  Ship  Research  and  Development  Center 


Bethesda,  Md.  2008L 

ATTN:  Code  56LI  1 

Commander 

Naval  Surface  Weapons  Center,  White  Oak 
Silver  Springs,  Md.  20910 

ATTN:  2L3  1 

2L0  1 

21*1  1 

Underwater  Explosions  Research  Division  1 


Naval  Ship  Research  and  Development  Center 
Norfolk  Naval  Shipyard 
Portsmouth,  Va.  23709 


Commander 

Naval  Surface  Weapons  Center,  Dahlgren  Laboratory 
Dahlgren,  Va.  22LL8 

ATTN:  MAL  I 

Common de r 

Naval  Facilities  Engineering  Command 
Department  of  the  Navy 
Hoffman  Building,  200  Stovall  Street 
Alexandria,  Va.  22332 

ATTN:  NFAC-OkB  1 

NFAC-03  1 

Chief  of  Naval  Research 
Navy  Department 
Arlington,  Va.  22217 
ATTN : Code  1*61* 
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Address 


No.  of 
Copies 


Navy  (Continued) 


Comman’ing  Officer 

Nuclear  Weapons  Training  Center 

Atlantic  Naval  Base 

Norfolk,  Va.  23511 

ATTN:  Nuclear  Warfare  Department 

Commander-In-Chief 
U.  S.  Atlantic  Fleet 
U.  S.  Naval  Base 
Norfolk,  Va.  23511 

Commander 

U.  S.  Naval  Oceanographic  Office 
Washington,  D.  C.  20373 
ATTN:  Library,  Code  1600 

Chief  of  Naval  Operations 
Navy  Department 
Washington,  D.  C 20350 
ATTN:  OP-75 

Director 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 
ATTN : Code  2627 


Special  Projects 
Navy  Department 
Washington.  D.  C. 
ATTN : CP-272 


20360 


Commander 

Naval  Ordnance  Systems  Command 
Navy  Department 
Washington.  D.  C.  20360 


Commander 

Naval  Ship  Engineering  Center 
Navy  Department 
Washington,  D.  C.  20360 
ATTN:  Code  6105 

Commanding  General 

Marine  Corps  Dev  and  Education  Command 
Quantico,  Va.  2213^ 

ATTN:  Director,  Development  Center 
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Navy  (Continued) 

Commandant  of  the  Marine  Corps 
Washington,  D.  C 20380 

ATTN:  Code  RD  2 


Air  Force 

Commander 

Space  and  Missile  Systems  Organization 
Norton  Air  Force  Base,  Calif.  92L09 

ATTN:  MNNH/CPT  J.  V.  Kaiser,  Jr.  2 

DEB  1 

Director,  USAF  Project  Rand 

Via:  U.  S.  Air  Force  Liaison  Officer 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  Calif.  90l06 

ATTN : Library  1 

Air  Force  Technical  Applications  Center/Tr  1 

Patrick  Air  Force  Base,  Fla.  32925 

Commander  1 

JSTPS/JLPW  (ATTN:  LCDR  M.  Barrata) 

Offutt  Air  Force  Base,  Nebr.  68113 


Commander 

Strategic  Air  Command 

Offutt  Air  Force  Base,  Nebr.  68113 

ATTN:  0AWS  1 

Commander 

Air  Force  Weapons  Laboratory,  AFSC 

Kirtland  Air  Force  Base,  N.  Mex.  87117 

ATTN:  DEV-G/M r.  R.  W.  Henny  1 

DEV-G/MAJ  J.  T.  Neal  1 

DbV-S/MAJ  James  M.  Warren  1 

DEV/CPT  G.  W.  Ullrich  1 

SUL,  Technical  Library  1 

DEV-S/Dr.  M.  A.  Plamondon  1 

Steve  Melzer  i 

DEV-F/Mr.  J.  L.  Bratton  1 


Air  Force  Institute  of  Technology 
AFIT  Bldg  6U0  Area  B 

Wright-Patterson  Air  Force  Base,  Ohio  15^33 

ATTN:  Technical  Library  1 


36 


Address 


No.  of 
Copies 


Air  Force  (Conv-inued) 


Commander  2 

Air  Force  Logistics  Command 
Wright-Patterson  Air  Force  Base,  Ohio  ^5^33 

Langley  Research  Center 
NASA,  Langley  Field 
Hampton,  Va.  23665 

ATTN : Mr . Philip  Donely  1 

Headquarters,  TAC/INTT  1 

Langley  Air  Force  Base,  Va.  23665 

Director  of  Civil  Engineering 
Headquarters,  Air  Force  Systems  Command 
Andrews  Air  Force  Base 
Washington,  D.  C.  20331 

ATTN : DEE  1 

Director  of  Civil  Engineering 
Headquarters,  USA? 

Washington,  D.  C.  20330 

ATTN:  AF/PREE  1 


APO  and  FPO 


Commander-in-Chief,  Pacific  1 

FPO  San  Francisco,  Calif.  9^129 

U.  S.  Documents  Officer  1 

Office  of  the  United  States  National 
Military  Representative-Shape 
APO  New  York,  N.  Y.  09055 


Domestic  Exchanges 


Northern  Arizona  University 
Box  15600 

Flagstaff,  Ariz.  86001 

ATTN:  Prof.  Sandor  Popovics  1 

Other  Government  Agencies 


Director 

Lawrence  Livermore  Laboratory 

P.  0.  Box  808,  Livermore,  Calif.  9^550 

ATTN:  K-Division 
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Other  Government  Agencies  (Continued) 

U.  S.  Department  of  the  Interior 
U.  S.  Geological  Survey 
3^5  Middlefield  Road 
Menlo  Park,  Calif.  9^025 
ATTN:  Mr.  Harold  W.  Olsen 

Director 

Lawrence  Livermore  Lahore * 

P.  0.  Box  808 
Livermore  Calif.  9^550 

Director 

Lawrence  Livermore  Laboratory 
Technical  Information  Division 
P.  0.  Box  808 
Livermore,  Calif.  9^550 
ATTN:  Technical  Library 

Ms.  Barbara  Germain 

Manager 

Albuquerque  Operations  Office 

U.  S.  Energy  Research  and  Development  Administration 
P.  0.  Box  5^00 
Albuquerque,  N.  Mex.  87115 

Commander,  Field  Command 
Defense  Nuclear  Agency 
Kirtland  Air  Force  Base,  N.  Mex.  87115 
ATTN:  FCTD,  Director 

FCPR 

Sandia  Laboratories 
P.  0.  Box  5800 
Albuquerque,  N.  Mex.  87115 
ATTN:  Technical  Library 

Commander 

Field  Command,  Defense  Nuclear  Agency 
Kirtland  Air  Force  Base,  II.  Mex.  87115 
ATTN : FCTD-T 

Sandia  Laboratories 

P.  0.  Box  5800 

Albuquerque,  N.  Mex.  87115 

ATTN:  Doc  Control  for  Mr.  L.  J.  Vortnan 
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Other  Government  Agencies  (Continued) 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663,  Los  Alamos,  N.  Mex.  875^ 

ATTN:  Document  Control  for  Reports  Library 


Sandia  Laboratories 

P.  0.  Box  5800,  Albuquerque , N.  Mex.  87115 
ATTN:  John  Kaiser 

Director 

Defense  Advanced  Research  Projects  Agency 
Architect  Building,  ll+00  Wilson  Blvd. 

Arlington,  Va.  22209 
ATTN : NTDO 

Administrator 

National  Aeronautics  and  Space  Adm. 

Washington,  D.  C.  205^6 

ATTN:  Marine  Transportation  Officer  (Code  BXC) 

National  Military  Command  Systems  Support  Center 
Pentagon  BE  685 
Washington,  D.  C.  20301 
ATTN:  Technical  Library 


Director  of  Defense  Research  and  Engineering 
Washington,  D.  C.  20301 
ATTN:  Technical  Library,  Room  3C-128 

Mr.  Frank  J.  Thomas 

Federal  Highway  Administration 
Department  of  Transportation 

Chief,  Structural  and  Applied  Mechanics  Division 
Washington,  D.  C.  20590 
ATTN:  Mr.  F.  J.  Tamanini 

Director 

Defense  Intelligence  Agency 
Department  of  Defense 
Washington,  D.  C.  20301 
ATTN:  DI-7D 

Director 

Defense  Civil  Preparedness  Agency 

Washington,  D.  C.  20301 

ATTN:  Mr.  George  Sisson  (RE-HV) 
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Other  Government  Agencies  (Continued) 

Director 

Defense  Nuclear  Agency 
Washington,  D.  C.  20305 
ATTN : SPSS 

Director 

Weapons  Systems  Evaluation  Group 
Washington,  D.  C.  20305 

U.  S.  Energy  Research  and  Development  Administration 
Technical  Information  Service 
Washington,  D.  C.  2051+5 

ATTN:  Chief,  Classified  Technical  Library 

Colleges  and  Universities 

University  of  Arizona 
College  of  Engineering 
Tucson,  Ariz.  85721 
ATTN:  Dr.  George  Howard 

University  of  Colorado 
School  of  Architecture 
Boulder,  Colo.  80301+ 

ATTN:  Prof.  G.  K.  Vetter 

University  of  Detroit 
Department  of  Civil  Engineering 
1+001  West  McNichols  Road 
Detroit,  Mich.  1*8221 
ATTN:  Prof.  Warren  J.  Baker 

George  Washington  University 

Nuclear  Defense  Design  Center 

School  of  Engineering  and  Applied  Science 

Washingtcn,  D.  C.  20052 

ATTN:  Prof  Raymond  R.  Fox 

George  Washington  University 

Nuclear  Defense  Design  Center 

School  of  Engineering  and  Applied  Science 

Washington.  D.  C.  20006 


1*0 


i — ~ -t-~-  • •»  ~ ^ 


No.  of 
Copies 


5 

1 


1 


1 


1 


1 


1 

1 


Address 


No.  of 
Copies 


Colleges  And  Universities  (Continued) 


University  of  Illinois 

Civil  Engineering  Building,  Urbana,  111.  6l801 

ATTN:  Metz  Reference  Room,  B106  1 

Prof.  W.  J.  Hall  1 

Prof.  A.  J.  Hendron,  Jr.  1 

Dr.  Nathan  M.  Newmark,  1211  1 

Iowa  State  University  of  Science  and  Technology 
231  Sweeney  Hall,  Ames,  Iowa  50010 

ATTN:  Prof.  Glenn  Murphy  1 


Lehigh  University 

Department  of  Civil  Engineering 

Bethlehem,  Pa.  18015 

ATTN:  Dr.  D.  A.  Van  Horn  1 

University  of  Massachusetts 
Department  of  Civil  Engineering 
Amherst,  Mass.  01002 

ATTN:  Dr.  M.  F.  White  1 

Massachusetts  Institute  of  Technology 

77  Masschusetts  Avenue,  Room  1-382,  Cambridge,  Mass.  02139 

ATTN:  Dr.  R.  V.  Whitman  1 

University  of  Michigan 

Department  of  Civil  Engineering 

30h  West  Engineering,  Ann  Arbor,  Mich.  48lOk 

ATTN:  Prof.  F.  E.  Richart,  Jr.  1 

University  of  Missouri  - Rolla 

Rock  Mechanics  and  Explosives  Research  Center 

Buehler  Building,  Rolla,  Mo.  65^01 

ATTN:  Dr.  George  B.  Clark,  Director  1 

University  of  New  Mexico  2 

Eric  H.  Wang  Civil  Engineering  Research  Facility 

P.  0.  Box  188,  University  Station,  Albuquerque,  N.  Mex.  87106 

Pennsylvania  State  University 
101  Eng.  A 

University  Park,  Pa,  16802 

ATTN:  Prof.  Richard  Kummer  1 
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Colleges  and  Universities  (Continued 

Pennsylvania  State  University 
Department  of  Architectural  Engineering 
University  Park,  Pa.  16802 
Prof.  G.  Albright,  Head 

Puriue  University 

School  of  Civil  Engineering 

Civil  Engineering  Building,  Lafayette,  Ind.  **7907 
Prof.  M.  B.  Scott 

Rensselaer  Polytechnic  Institute 
Mason  House,  Troy,  N.  Y.  12180 

ATTN:  Dr.  Clayton  Oliver  Dohrenwend,  Security  Officer 

Rice  University 

Department  of  Civil  Engineering 
P.  0.  Box  1892,  Houston,  Tex.  77001 
ATTN:  Prof.  A.  S.  Veletsos 

San  Jose  State  College 
Department  of  Civil  Engineering 
San  Jose,  Calif.  9511** 

ATTN:  Dr.  Robert  G.  Spicher 

University  of  Texas  at  Austin 
Balccnes  Research  Center 
Route  U,  Box  l89»  Austin  Tex.  78757 
ATTN:  Dr.  J.  Neils  Thompson 

Utah  State  University 

Department  of  Mechanical  Engineering 

Logan,  Utah  8** 321 

ATTN:  Prof.  R.  K.  Watkins 

University  of  Washington 

Department  of  Civil  Engineering  FX-10,  Seattle,  Wash. 
ATTN : Chairman 

Worcester  Polytechnic  Institute 
Department  of  Civil  Engineering 
Worcester,  Mass.  01609 
ATTN:  Dr.  Carl  Koontz 
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Corporations 


Aerospace  Corporation 

P.  0.  Box  92957*  Los  Angeles,  Calif.  90009 

ATTN:  Dr.  Prem  N.  Mathur  1 

Agbabian  Associates  2 

Fin  nppr'i  nrr  Consultants 

250  North  Nash  Street,  El  Segundo,  Calif.  9021+5 
Applied  Theory  Incorporated 

1010  Westwood  Boulevard,  Los  Angeles,  Calif.  90021+ 

ATTN:  Dr.  John  G.  Trulio  1 

AVC0  Corporation 

Research  and  Advanced  Development  Division 
201  Lowell  Street,  Wilmington,  Mass.  01887 


ATTN:  Mr.  R.  E.  Cooper  1 

Bell  Telephone  Laboratories 
Outside  Structures  Department 
Whippany,  N.  J.  07981 

ATTN:  Mr.  John  Foss  1 

The  Boeing  Company 

P.  0.  Box  3707,  Seattle  Wash.  98121+ 

ATTN:  Technical  Library  1 

Battelle  Columbus  Laboratories 

505  King  Avenue,  Columbus,  Ohio  1+3201 

ATTN:  Dr.  P.  N.  Lamori  1 

Energy  Incorporated  1 

P.  0.  Box  736,  Idaho  Falls,  Idaho  831*01 

General  Research  Corporation  1 

Document  Control  Supervisor 


Westgate  Research  Park,  McLean,  Va.  22101 

General  Electric  Company,  Tempo 

8l6  State  Street,  Santa  Barbara,  Calif.  93102 

ATTN:  Mr.  Warren  Chan  (DASIAC)  1 

General  Electric  Company 

Space  and  RESD  Division 

P.  0.  Box  8555,  Philadelphia,  Pa.  19101 

ATTN:  Larry  Chasen,  Manager-Libr . 1 
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Corporations  (Continued) 

General  Research  Corporation 
P.  0.  Box  3587 

Santa  Barbara,  Calif.  93105 

ATTN:  Mr.  Ben j amine  Alexander  1 

Dr.  John  S.  Rinehart  1 

Hyperdynamic s 

P.  0.  Box  392,  Santa  Fe,  N.  Mex.  87501 
I IT  Research  Institute 

10  West  35th  Street,  Chicago,  111.  60616 

ATTN:  Dr.  T.  Schiffman  1 

Martin  Marietta  Aerospace,  Orlando  Division 
P.  0.  Box  5837,  Orlando,  Fla.  32805 

ATTN:  Gerbert  E.  McQuaig  (MP  18)  1 

Physics  International  Company 

2700  Merced  Street,  San  Leandro,  Calif.  9^577 

ATTN:  DOC  Control  for  Dr.  Charles  Godfrey  1 

DOC  Control  for  Mr.  Fred  M.  Sauer  1 

R&D  Associates 

P.  0.  Box  9695 > Marina  Del  Rey,  Calif.  90291 

ATTN:  Dr.  Harold  L.  Brode  1 

Dr.  H.  F.  Cooper,  Jr.  1 

Dr.  Bruce  Hartenbaum  1 

Mr.  J.  G.  Levis  1 

The  Rand  Corporation  1 

1700  Main  Street 

Santa  Monica,  Calif.  90406 

Southwest  Research  Institute 

8500  Culebra  Road,  San  Antonio,  Tex.  78228 

ATTN:  Dr.  Robert  C.  Dehart  1 

Systems,  Science  and  Software 

P.  0.  Box  1620,  Le  Jolla,  Calif.  92307 

ATTN:  Dr.  Donald  R.  Grine  1 

Dr.  T.  R.  Blake  1 

Electromechanical  Systems  of  New  Mexico,  Inc. 

P 0.  Box  11730,  Albuquerque,  N.  Mex.  87112 

A^TN:  R.  A.  Shunk  1 
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Corporations  (Continued) 

Stanford  Research  Institute 
333  Raven swood  Avenue 
Menlo  Park,  Calif.  91*025 

ATTN:  Dr.  G.  Abrahamson  1 

Terra  Tek,  Inc. 

University  Research  Park 

1*20  Wakara  Way 

Salt  Lake  City,  Utah  81*108 
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In  accordance  with  Bt  70-2-3,  paragraph  6c(l)(b), 
dated  15  February  1973,  * facsimile  catalog  card 
la  Library  of  Ccograaa  fonat  la  reproduced  below. 


Drake,  Janes  L 

A method  for  designing  deep  underground  structures  sub- 
jected to  dynamic  loads,  by  Janes  L.  Drake  and  James  R. 
Britt.  Vicksburg,  U.  S.  Army  Engineer  Waterways  Experi- 
ment Station,  1976. 

45  p.  illus.  27  cm.  (U.  S.  Waterways  Experiment 
Station.  Technical  report  N-76-9) 

Prepared  for  Defense  Nuclear  Agency,  Washington,  D.  C., 
and  Office,  Chief  of  Engineers,  U.  S.  Army,  Washington, 
b.  C.,  under  DNA  Subtask  J34CAXSX311,  "Underground  Struc- 
tures Studies,"  and  OCE  Project  4A762719AT40/A1/017, 
"Stability  of  Deep  Underground  Structures  In  Rock." 

References:  p.  27. 

1.  Dynamic  loads.  2.  Underground  structures. 

I.  Britt,  James  R. , joint  author.  II.  Defense  Nuclear 
Agency.  III.  U.  S.  Army.  Corps  of  Engineers.  (Series: 
U.  S.  Waterways  Experiment  Station,  Vicksburg,  Miss. 
Technical  report  N-76-9) 

TA7.W34  no  N-76-9 


